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We study the localization properties of the anomalous diflusion phase x ^ t'^ with < /i < 
1 which exists both in the Sinai diffusion at small bias, and in the related directed trap model 
presenting a large distribution of trapping time p(t) ~ l/r^^*^. Our starting point is the Real 
Space Renormalization method in which the whole thermal packet is considered to be in the same 
renormalized valley at large time : this assumption is asymptotically exact only in the limit of 
vanishing bias and corresponds to the Golosov localization. For finite ^, we thus generalize 
the usual RSRG method to allow for the spreading of the thermal packet over many renormalized 
valleys. Our construction allows to compute exact series expansions in ^ for all observables : to 
compute observables at order /x", it is sufficient to consider in each sample a spreading of the thermal 
packet onto at most (1 + n) traps. So our approach provides a description of the structure of the 
thermal packet sample by sample, and a full statistical characterization of the important traps at 
a given order in /i. For the directed trap model, we show explicitly up to order /i^ how to recover 
the exact expressions for the diffusion front, the thermal width, and the localization parameter ¥2- 
We then use our method to derive new exact results for the localization parameters Yk for arbitrary 
k, the correlation function of two particles, and the generating function of thermal cumulants. We 
then explain how these results apply to the Sinai diffusion with bias, by deriving the quantitative 
mapping between the large-scale renormalized descriptions of the two models. We finally study the 
internal structure of the effective 'traps' for the Sinai model via path-integral methods. 
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I. INTRODUCTION 

The motivations to study the Sinai model Q] have two origins. On one hand, the Sinai model represents a 'toy' 
disordered system, in which many properties that exist in more complex systems can be studied exactly, such as 
for instance aging behaviors ^ , persistence exponents |^ , the decoupling of the dynamics into fast degrees of 
freedom which rapidly reach local equilibrium and a slow non-equilibrium part governed by metastable states 
and some chaos and rejuvenation effects Q. On the other hand, the Sinai model directly appears in various specific 
systems, ranging from the dynamics of domain walls in the random field Ising chain (7[ ||] to the unzipping transition 
in DNA ^ . It is thus interesting to obtain exact detailed informations for various observables in the Sinai model. 

One of the most important property of the symmetric Sinai diffusion is the following localization phenomenon 
discovered by Golosov jl^ : all the thermal trajectories starting from the same initial condition in the same sample 
remain within a finite distance of each other even in the limit of infinite time! In particular, in a given sample, for 
a given initial condition, the rescaled position X — -f0^ is deterministic, and it is only after averaging over the 
samples that X is distributed with the Kesten distribution j|, |o[ 0. The physical picture is that the particle is at 
time t near the bottom of the deepest valley it has been able to reach. This is why the Real Space Renormalization 
Group method, first introduced in the field of random quantum spin chains [p^ [l^ , is so well suited to study the 
symmetric Sinai diffusion Q. Recently we have studied in more details the localization properties, by computing 
the infinite time limit of the localization parameters, which represent the disorder-averages of the probabilities that 
k independent particles in the same sample starting from the same initial condition are at the same place at time t 
and of the correlation function C{l,t), which represents the disorder- average of the probability that two independent 
particles in the same sample starting from the same initial condition are at a distance I from each other at time t. We 
have moreover shown ||^ that the the infinite time limit of the localization parameters and of the correlation function 
exactly coincide with the corresponding equilibrium observables in a Brownian potential in the thermodynamic limit. 

A natural question is thus : do some of these localization properties survive in the presence of a small bias ? 



A. Sinai model with bias 



The Sinai model in the presence of a constant bias Fq > can be studied in a continuum version via the following 
Langevin equation 

^=Fo~U'ixit))+riit) (1) 

where r]{t) is the usual thermal noise 

<r]{t)r]{t') >=2TS{t-t') (2) 
and where U{x) is a Brownian random potential representing the disordered landscape 



{U{x) - U{y) f = 2a5{x - y) (3) 

Equivalently, the model may be defined by the Fokker-Planck equation in a given sample {[/(x)} for the probability 
distribution P{x, t\xQ, 0) 

dtP{x, i|xo, 0) = [r9, + U'{x) - Fo] P{x, t|xo, 0) (4) 

In the biased case Fq > 0, the diffusion becomes transient, and there are dynamic phase transitions |l^ as 

fo grows, in terms of the dimensionless parameter 

= (5) 

a 

For < /i < 1, the mean position of the particle presents the anomalous behavior ||lj, |l^, |l^ 



< x{t) > cx t'" (6) 

t — ^oo 



whereas for fi > I, there is a finite velocity < x{t) > ^ V{fi)t. For the anomalous diffusion phase, the exact diffusion 
front is given in terms of Levy stable distributions ||lj, |l5[ 17, |l^: we refer the reader to the Appendix ^ for the 
definition and properties of these Levy fronts. 
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B. Directed trap model 

It has been suggested in that at large time, the physics of the Sinai model with bias is actually equivalent to a 
simple directed trap model, defined by the Master equation 

dPtjn) ^ _Pt(n)_ ^ Pt{n-1) 
dt T„ r„_i 

with the initial condition Pt^oin) = Snfi, and where the trapping times are independent random variables distributed 
with a law presenting the algebraic decay 

« i (8) 

The anomalous diffusion phase < < 1 then corresponds to the phase where the mean trapping time < t >— 
J dTTq{T) is infinite. The corresponding diffusion front is also a Levy diffusion front (see Appendix ^) as for the 
biased Sinai diffusion discussed above. For simplicity in the article, we choose the normalization of the algebraic tail 
to be 

Qir) ^ (9) 

It is clear that this choice simply amounts to a rescaling of t. 

The presence of the large algebraic decay in the effective trapping time distribution (^) for the biased Sinai diffusion 
may be understood from the Real Space Renormalization approach in relation with the distribution of the barriers 
against the drift in the renormalized landscape at scale T ^ 

Pr{F) = e{F - r)2^e-2*(-^~^^ (10) 

where 2S = The trapping time r ~ e^^ is then distributed with the power law (^) with the correspondence 
^ = 2ST. 



C. Previous results for the localization in the directed trap model 

For the directed trap model, the existing results on the extension of the thermal packet are twofold. On one hand, 
the thermal width has been exactly computed in [19] (equation (26)) 



+00 +00 / . s 3 



<An2(t)>^^n2F,(n)~[^nP,(n)P = -i-(^) 7(^)^2^ (n) 

where the integral /(/i) of Equation (26) in jl^ can be rewritten after a change of variables as 

IM I'd. (1 + ^)^"^-^)^" (12) 
l(^i) - ''%2^+2 + 2cos7r/.zA<+i + l ^ ' 

The result ( pi] ) shows that the the thermal packet is spread over a length of order t^^. 

On the other hand, the infinite-time limit of the localization parameter for fc = 2 has been exactly in |2^ : their 
result (24) may be rewritten after a deformation of the contour in the complex plane as 

+ °° f^'^ df) piSfi _ piO 

n=0 -^-^ 

This expression shows that Y2 is finite in the full phase < fi < 1 and vanishes in the limit fi — 1. How can this 
property coexist with the result (|l^) for the thermal width ? The numerical simulations of pO[ | show that for a single 
sample at fixed t, the probability distribution Pt(n) is made out of a few sharp peaks that have a finite weight but 
that are at a distance of order f^. This explains why at the same time, there is a finite probability to find two particles 
at the same site even at infinite time, even if the thermal width diverges as at large time. 
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D. Goal and results 



The aim of this article is to provide a probabilistic description, sample by sample, of the localization properties of 
the directed trap model and of the Sinai diffusion with bias in the anomalous diffusion phase < /i < 1. We will 
need to generalize the usual Real Space Renormalization group method to allow for the spreading of the thermal 
packet over many renormalized valleys. Indeed, in the usual RSRG method, the whole thermal packet is considered 
to be in the same renormalized valley at large time : this assumption is asymptotically exact in the symmetric Sinai 
model and actually corresponds to the Golosov localization ||^, |l^ discussed above; it is also valid for the biased case 
but only in the double limit of vanishing bias /i — > and large time with the fixed parameter 7 = jjT \nt [|j . We will 
thus define explicit rules for the RSRG approach with multiple valleys occupancies and show that our construction 
allows to compute exact expansions in /i for all observables. 



1. Summary of results for the direeted trap model 



For the directed trap model, we explicitly show how to recover in a unified framework the expansions up to order 
of the exact results for the observables discussed above : 

• Expansion in /i of the Levy diffusion front for the rescaled variable X = (see Appendix ^) 



g{X) =e-^+A*7B(^-l)e 



-X 



12 



^ + —\ +X\—-?.^\+X'\^- — 



12 



? 9 

.2 / Te 7^ 



12 



-X 



(14) 



Expansion in fi of the thermal width ( from ( [TTI ) and ([T 



A(m) = lim 



< An2(t) > 



t2, =M21n2)+M^[-- 



2 In 2(ln 2 - 2 + 27^)] + 0{fi^) 



(15) 



Expansion in fi of the localization parameter Y2 (from dl 



Y2{p) = 1 - /i(21n2) + /i2(41n2 - — ) + O(Ai') 
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(16) 



These comparison with exact results show that our generalized RSRG procedure is exact order by order in /j, : to 
compute observables at order it is sufficient to consider a spreading of the thermal packet onto at most (1 + n) 
traps. So our description provide a description of the structure of the thermal packet sample by sample, and a full 
statistical characterization of the important traps at a given order in fi. 

We then use our procedure to derive new exact results. We obtain the expansion in fi of the localization parameter 
Yk for arbitrary k up to order /i^ 



Ykifi) - l + M 



dv , 



-kv 



/■+°° rlv /■+°° din 

+^2 / _ / ^[pl{^y^^)+pl(^yj^y) + 2pl{v,w) + l-e-'^'^ -2{l-e-^f-[l-e-^f] + 0[ii^) (17) 
Jo V w 



where the functions p2 and are defined in ( \8i 

We obtain that the correlation function of two particles averages over the disorder reads 



+00 +00 ^ y , 

C{l,t) = V VF(n,t|0,0)F(TO,<|0,0),5,.|„_„.| ~ Y2{fi)6i^o + -cj - 



n—0 7n—0 



(18) 



where the weight of the delta pic at the origin corresponds as it should to the localization parameter Y2 (p^), whereas 
the second part presents a scaling form of the variable X ~ We obtain the following expansion for the scaling 
function C,, 



Cm (A) 



-A 



^(2 In 2) + 



— -ln2(4 + ln2 + 7£;) + A ( -y + In 2(ln 2 + 7^) 



(19) 
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We also consider the generating function of rescaled thermal cumulants Cfc(/i) 



+00 



Zfj,{s) = In < e^""^ > = X! 



k=l 



kl 



■c-k(p) 



The first one simply represents the mean value which can be obtained from the diffusion front (|lj 



(20) 



(21) 



The second one C2{n) represents the thermal width A(/i) (|15[). We obtain the expansion at first order in /i of the 
generating function 



-s + fj, 



/ dFe-^ / - In [e-^ + (1 - e-'')e-^^] + / - In [e-"e^^ + (1 - g-")] 
Jo Uo ^ Jl V 



Oil?) (22) 



The series expansion in s then yields all thermal cumulants at first order in /i. In particular, the first terms beyond 
the mean value Ci {fi) and the thermal width C2 (/i) read 



C3(/i) = lim 

t — *oo 



< rr' > -3 < >< n > +2 < n >^ 



/i6(21n3- 3 In 2) + 0(/i2 



C4(/i) 



lim 

t — ^oo 



< n'^ > -A < rfi >< n > -3 < n2 >2 +12 < >< n >2 -6 < n 



t2M 



/Lt24(191n2 - 12 In 3) + ©(/Lt^) 



(23) 



(24) 



^. Summary of results for the biased Sinai model 

We will derive an exact quantitative mapping between the renormalized descriptions of the trap model and the 
biased Sinai diffusion with bias. As a consequence, in the whole anomalous diffusion phase < /i < 1, all properties 
of the directed trap models that concern the rescaled quantity X — ^ are exactly the same for the Sinai model with 
bias in terms of the rescaled quantity 



This relation was already guessed in ||lj] for the special case of the averaged diffusion fronts of the two models (see 
Appendix 0). In particular, the thermal width of the Sinai model reads from the exact result ( ITTl) of ||l^ 

<AxHt)> _{a^(3'y'TH^^) fsinn^,Y ^^^^ ^^g) 

(27) 



t2M cr2^4 Y{2n) V TT/i 



2fl3\2M 



a2/34 



^^ + h^ + 21n2(ln2-2-27s)]4 + 0(-) 



and more generally, all thermal cumulants can be obtained from the results of the trap model ( p^ via the correspon- 
dence (ESh. 



For the localization parameters, the result yj*''"? represents for the biased Sinai model the probability to find k 
independent particles at a finite distance of each other in the limit of infinite time. These particles are then distributed 
with the Boltzmann distribution in a infinitely deep biased Brownian valley, leading to 

Ysinai yrtrapyrvalley (28) 



where Y^"" "^^j computed in ( 209 ), is the localization parameter for k particles at equilibrium in a infinitely deep biased 



Brownian valley. 

For the two-point correlation function, we obtain for the biased Sinai model the two-scaling form 

Cs^na^{l,t) = Y^"-"-^ C,aUey{l) + ^2 (^) [ta^ ^^f^ ^ " r2 (^) [tCT^ (}^f^ ^^^^ 

where the first part represents the case where the two particles are at a finite distance of each other at infinite time, 



in which case their correlation Cyaiieyil) is given by (211). The second part, corresponding to the cases where the 
two particles are in different renormalized valleys at infinite time, is exactly given by the scaling function ( [l9| ) 
describing the long-range behavior in the trap model. 
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3. Organization of the paper 



We first study the directed trap model: the Section || presents the usual RSRG which yields all observables in the 
limit /I — !■ 0; in Section III we explain the origin of the spreading of the thermal packet at first order in n and compute 
observables at this order ; in Section IV we study the second order fj,^ ; in Section ^ we explain the structure of the 
set of important traps at any given order /i". 

We then turn to the biased Sinai model: in Section Vl, we derive the quantitative mapping between the large-scale 



renormalized descriptions of the two models (the biased Sinai model and the directed trap model); in Section VII, we 
moreover characterize the internal structure of the 'traps' in t he bi ased Sinai model by computing various statistical 
properties of infinitely biased Brownian valleys. The Section VIIl contains a discussion of the universality. Finally, 
the Section IX contains the conclusion, and some more technical details are given in the Appendices. 



II. DIRECTED TRAP MODEL IN THE LIMIT /i ^ 

The Real Space Renormalization procedure for the Sinai model ||^ can be reformulated for the directed trap model 
as follows. At time t, all traps with trapping time Ti < t are decimated and replaced by a "flat landscape" to produce 
the renormalized landscape at time t. We stress here that, contrary to the symmetric Sinai diffusion, the remaining 
traps are just some of the initial traps, and that their trapping time have not been renormalized by the decimation of 
the small ones. This non-renormalization of the trapping times actually corresponds for the biased Sinai landscape to 
the fact that barriers against the bias converge without rescaling to a fixed distribution . The usual RSRG picture 
for the dynamics is now very simple : the particle starting at t = in the n — trap will be at time t in the first trap 
of the renormalized landscape, that is in the first trap having a trapping time bigger than t. We will call this trap 
the Main Trap M. In the usual RSRG approach, all thermal trajectories are all in the same trap M. In particular, the 
probability distribution in a given sample is a delta 

pW(n)=<5„,.„ (30) 

and the localization is total : there are no thermal fluctuations 

[An2(t)](0) = (31) 

and more generally all thermal cumulants beyond the first one vanish : the generating function of thermal cumulants 
( pO| ) simply reads 

Z(°)(.) -~s^ = -s (32) 
The two-particle correlation function is a delta 

+00 4-oc 

C^°Hl,t) ^J2T. P^°\n)P^°\m)Si,^„_^^ = Si,o (33) 

n— m— 

and the localization parameters have their maximal value 

Y^"\t)^l (34) 

The corresponding averaged diffusion front is thus simply given by the distribution of the position n = um of the 
main trap 



P(")(n) 



+ 00 



1 - / drq(T) 



+ 00 



drq(T) (35) 



where the first part [..]" represents the probability that the first n traps have a trapping time < t, and where the 
last part represents the probability that the n**^ trap has a trapping time > t. So the scaling function g describing 
the averaged diffusion front at large time 

1 / n 
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FIG. 1: Hierarchical structure of the important traps for a particle starting at the origin. The dashed line separates the 
"small" traps (that have a trapping time smaller than t) from the "big" traps (that have a trapping time bigger than t). The 
first big trap called M is occupied with a weight of order 0{fjP). The next big trap Li and the biggest small trap Si before M 
are occupied with weights of order 0{fi). The third big trap L2, the biggest small trap I2 between M and Li, and the second 
biggest small trap 52 before M are occupied with weights of order 0(/i^). 



is given at this order by a simple exponential 

5(0)(X) = e-^ (37) 

which indeed coincides with the limit /i — s- of the exact Levy front (see Appendix ^). 

So the approximation where all particles of the same thermal packet are in the same trap is correct only in the 
limit of vanishing fi. For finite fj,, we will have to allow for a possible dispersion of the thermal packet. In fact in the 
limit II we have considered that the distribution of the trapping time was infinitely broad in the following sense 
: all traps with < t were such that ^ ^ 0, whereas all traps with Ti > t were such that ^ ~ +00. For finite n, we 
have to take into account that these ratios are not really zero or infinite. We will do it order by order in fi. 



III. DIRECTED TRAP MODEL AT FIRST ORDER IN fj, 
A. Origins of the dispersion of the thermal packet at order /i 

At first order in /i, we need to consider two effects (see Figure |l|) : 

• the main trap M defined above has a trapping time tm which is not infinite. There is a small probability 

_ t 

(1 — e ) that the particle has already escaped from this main trap Af at time t, to jump into the next 
renormalized trap that we will call Li (for Large trap number 1), which is defined as the second trap satisfying 

n > t. 

• the biggest trap before the main trap, that we will call Si (for Small trap number 1), has a trapping time rg^ < t 

_ t 

which is not zero and thus there is a small probability e ''^1 that the particle is still trapped in Si at time t. 
We now describe the statistical properties of these two effects. 
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B. Statistical properties of the trap Li 

The joint distribution of the trapping time tm and of the position n — um of the main trap M and of the position 
U]^ of the next renormahzed trap Li reads 



dTq{T){m) 





P + OO 


n 


P + CO 


riL-n-l 


T^M,Li {n, TiL] Tm) = 0(t < TM)0(n < ul) 


1 - drqir) 




1- drqir) 


Jt 



with the scahng function 



2^-t^J\/,Li — 77"' — —iTm 



Dm.lAX,Xl;tm) = Oit < TM)0iO < X < Xl] 



(39) 



In particular, the distribution of the trapping time tm is obtained, as it should, by simply normalizing the original 
distribution (/(r) on the interval [t, +oo[ 



qt{TM) = fdX f dXLDM,LAX,XL;TM) = 0{tm - 



(40) 



The probability ttli = ^1 — e ^ to have already escaped from the main trap at time t and to be thus already ir 
the trap Li reads after averaging over the disorder 



dTqtir) 
1 dv 



1 - e- 



) =Ai /" ^^^"^(1-6-") 



1 dv 



■lnt;(l-e-") + 0(/i3) 



so it is of order fi. 

At this level of approximation, the diffusion front for a given sample is made out of two delta distributions 



(41) 
(42) 

(43) 



C. Statistical properties of the trap Si 

The trap Si has been defined as the biggest trap before the main trap M. The joint distribution of the position n 
of the main trap, the position ng and the trapping time Tg of the trap 5*1 read 



1^Si,M(.'^S,n;Ts) =9{t>Ts) 



1 - / drqir) 



4-00 



q{r,) / drqir) (44) 



LDs,^M{Xs = ^,X^^^;rs) (45) 



t— »oo t-»oo 

where the scaling function reads 

D s,.MiXs,X-rs)^ Oit > rs > l)eiX > Xs>Q)^\ — ] e~^y~) (46) 

TS \TS J 

We note that here, there are correlations between the trapping time and the positions, contrary to the decoupled 
measure ( ^9|) concerning the trap Li. The joint distribution of the positions alone reads 

SiXs,X)^JdrsDs,MXs,X;rs) = &Z|5>0)^_^ ^^^^ 

i.e. X is distributed with = (^7|), and Xs is uniformly distributed on the interval [0,X]. The distribution 
of the trapping time alone reads 



r 



+00 P+OO r^~'^ 

(ts)= / dXDs„MiXs,X;rs)^9it>rs)fi^ (48) 



The probability ns^ = e to be still in the trap 5*1 at time t reads after averaging over the disorder 

' dv 



dTsr{Ts)e = n 



/i V 



Ji V 



so it is of order ^. 

At this level of approximation, the probability distribution reads 

^Sim'^^(") = e~^<5„,„si + (l - e"^) (5„,„„ 
We now use the statistical properties of the traps Li and to compute various observables at order /i. 



(49) 
(50) 



(51) 



D. Diffusion front at order /i 

The correction due to the trap Li to the diffusion front in a given sample ( ^3|) with respect to one delta function 
at the zero-th order (pO|) reads 



_ d(o)+(i). 



(52) 



The average over the samples, that is over the positions {nM,nL) and over the trapping time t with the measure 
(B9) yields the correction to the scaling function (p6h 



dTM / dX / dXLDMM {X, Xl; tm) 1 - e ^ (,5(r - X^) - 5(y - X)) 



iX-l)n / dw^-i(l-e-") 



'(X-1) 



e-") +^? C— In 1,(1 - e-'') + 0(^.3) 
. Jo Jo ^ 



(53) 
(54) 
(55) 



Similarly, the correction due to the trap Si to the diffusion front in a given sample (|51|) with respect to one delta 
function at the zero-th order ( |30| ) reads 



After averaging over the samples with the measure ( ^6| ) , we obtain the correction to the scaling function ( |3^ ) 
5s?m(^) ^ I dTs I dX I dXsDs,M{Xs, X; Ts)e-^[S{Y - Xs) - <5(r - X)] 



/ ^^!!e-"(l-ri;^)e 
Ji w 



/■+°° r~^°° dv 

/ie-^(l-X) / ^e-''+/i'e-^(^'-2X) / _ Int-g-^ + 0(^3) 
Jl V Ji v 



Adding these two contributions to the zeroth order front (|37|), we finally get 

a'TJ'' ^ 9^'\X) + .9« ^ {X) + gi%{X) = e-^ + e-^{X - 1)^7^. + O(m^) 
which coincides with the expansion at order jj, of the exact Levy front (|lj). 



(56) 

(57) 
(58) 
(59) 

(60) 
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E. Thermal width at order fi 



For a given sample, the contribution of the trap Li to the thermal width reads (43) 

[< An^{t) >]''^^^^ =<7i2>-<n>2== e"^(l - e^^){nL ~ nuf (61) 

Averaging over the disorder, that is over the positions and over the trapping time r with the measure ( p^ ) yields : 

(1) 



< An2(i) > 



ML, 



2^l I dvv^'-^e-''{l-e-'') 
Jo 



(62) 



^ 2^ C — e-"(l - e-") + 2^? C — ln^;e-''(l - e"") + 0{^?) (63) 

Jq V Jo V 

Similarly, the contribution of the trap Si ( |5l| ) averaged over the samples with the measure (^6|) reads 

-1 (1) 



M 



< An2(t) > 



J SiM 



e-^(l_e-^)^^^^^^^-2;U dvv-'-^'^e-^l-e-^) (64) 



1 

+ 00 



= 2/i/ — e-''(l-e-'')-6^M _ in^;e-"(l _ e"") + 0(^3) (gg) 
Ji V Ji V 

Adding the two contributions finally yields 

[A(A^)]|oL = 2^ / - e"'') - M(21n2) (66) 

in agreement with the exact result ( |T5| ) of [p^ . 



F. LocaHzation parameters at order ^ 



For a given sample, the contribution of the trap Li (^3|) to the localization parameter Yk representing the probability 
to find k independent particles in the same trap at time t reads 



1 — e 



(67) 



So after averaging over the samples, that is over the trapping time r (^0|), the correction to the zero-th order ( |34| ) 
due to the trap Li reads 



ml>^^{t) ^ [Y,]^!it^ \i) - Yt \t) - M / dvv^~' [e-^^ + (1 - e-^'f 1 



(68) 



Similarly, the correction to the zero-th order ( p^ ) due to the trap 5*1 reads after averaging over the samples, that 
is over the trapping time Tg (Q) 



[y,]« ,(i) ^ mfj}'\t) ~ Yl^'\t) = M / dvv-^-^[e-^^ + (1 - e-^-f 1] 



Adding these two contributions, we finally get at first order in 

C-\-OQ 



y(0)+(i) ^ ^(0) ^ [^^j W ^ + ^ 1 + ^ p ^[e-'^^ + (1 - - 1] + O(M^) 



For the special case k — 2, the result 



yW+W = l_(21n2)/i + 0(M') 
is again in agreement with the expansion (O) of the exact result pOj . The other first values of k yield 



(o)+(i) 



Y. 



= 1 - (31n2)^ + 0(/i2) 



Y 



(o)+(i) _ 



32. 



= 1 - (21n— )/i + 0(^2) 



(69) 
(70) 

(71) 

(72) 
(73) 
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G. Correlation function of two particles at order fj, 

We now consider the correlation function of two particles (p^). For a given sample, the contribution of the trap Li 
(p3|) reads 



-oo +00 



■m\ 



(74) 



(75) 



After averaging over the disorder with the measure (|3^), the correction with respect to the zero-th order of the 
correlation function (p3|) reads 



1 



It presents the form ( p^ ) : the weight of the 5 part has been obtained in (^) and the scaling function reads 

Jo 

2^ /' ^e-(l _ e-") + 2a*' /' — ln^;e-''(l - e"") + 0(/? 
. Jo V Jq V 



(76) 

(77) 
(78) 



Similarly, the contribution of the trap Si leads after averaging over the samples with the measure (|4q ) 



4'm(M) ^ 4°M - C(0)(/,0 = [(e-^)' + (1 - e-^y l]<5,o + 2^-^ (l - e"^) 5,(„„_„,,) 



iW x.^^o.. I ^e-''(l-e-")le-^^^ 



It presents the form ( p^ ) : the weight of the S part has been obtained in ( |69| ) and the scaling function reads 



+00 



dv 



2/1 



dv /■+°° (7?) 

— e-"(l - e-") - 2/.2A / — lnz;e-"(l - e~n 
V .h V 



(79) 

(80) 
(81) 



Adding these two contributions, we finally get at first order in fi the following correlation function 



[1 - (21n2)/i + 0(/i2)] 5,^0 + Le-w [(2 In 2)/^ + 0{ii^ 



(82) 



H. Generating function of thermal cumulants at first order in 



The correction to the generating function ( |2C| ) due to the trap Li ( |43| ) with respect to the zero-th order ( |32| ) reads 
with the measure 



[Z^(.)]«^^ ^ [Z,{s)]^^}+^''^ - Zf\s) = In + (1 - e-^)e-(^--^) 



+00 



dw^-i / dFe-^ In [e-" + (1 - e-")e 



(83) 



Similarly, the correction due to the trap 5*1 reads with the measure 



[Z,(.)]W , ^ [Z,{s)]fJ''^ - Zl^\s) = In 
-°° dv '•+°° 



— / dre-^" In [e-"e^^ + (1 - e-") 



(84) 



12 



The total correction at order /i thus reads 



dYe- 



f — In [e-" + (1 - e-'')e-^^] + — In [e-"e^^ + (1 - e"")] 



(85) 



We may now perform a series expansion in s and evaluate the integrals to obtain the generating function of all thermal 
cumulants at first order in /i 



total 



-s-iE + s^ln2 - s3(21n3 - 31n2) + s'*(191n2 - 12 In 3) + 0(5^) 



(86) 



leading to the results (pj). 



IV. DIRECTED TRAP MODEL AT ORDER 

A. Dispersion of the thermal packet at order /i^ 

To compute observables at order /i^, we now have to consider the possible dispersions of the thermal packet over 
three traps. Denoting by ti and T2 the first two trapping times, the occupation probabilities of the three ordered sites 
are given by 

Pi{t]Ti) ^e^~ 
P2{t;Ti,T2) = e -2 -e " 

T2 - Tl V 

t ^ 

P3{t;n,T2) =1-^ (87) 

T2 - Tl 

In the following, we will also use the notations v = —, w = — 

Pi{v) ^e-- 

P2{v,w) =^^(e— -e-) 

V — w 

P3{v,w) = 1 - pi{v) - p2{v,w) (88) 

To simplify computations later, it will be convenient to use in intermediate calculations the two following obvious 
properties : the occupation probability of the third site is a symmetric function of {v, w) 

P3{v,w) = P3{w,v) (89) 

and the three occupation probabilities satisfy the normalization 

Piiv) +p2iv,w) + P3{v,w) = 1 (90) 

The enumeration of the various possibilities for the three traps is as follows (see Figure if) . 



Configurations [M, L\ , L2 



The tree traps are the main trap M, the next renormalized trap Li introduced in (IIIB) and the second next 
renormalized trap that we call L2. The joint distribution of the rescaled positions and trapping times read 



Tm,LulAX,Xml,,Xl,;tm,tl,) = 9(1 < TM)0{t < tl,)9{0 < X < Xi < X2) 



Til \TLi 



At this level of approximation, the diffusion front is made out of three delta pics 

,(0) + (l) + (2)/ X _ 



ML1L2 



(n) = pi{t]TM)5 + P2it;TM,TL,)S + P3it;TM,TLi)5 



(91) 



(92) 



where the weights of the three traps are given by dS 
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2. 



Configurations {Ad, I2, Li) 



The tree traps are the mam trap M, the next renormahzed trap Li mtroduced in (IIIB) and in between the 
intermediate trap that we call /2, defined as the biggest trap in the decimated region between M and Li. The joint 
distribution of the rescaled positions and trapping times read 



TMJ,,L^ {X, Xi,Xl; tm,ti) = 9{tm > t > n > 1)9{Xl > Xi > X > 0) 



(1 

Tl \TI 



The corresponding diffusion front reads 



(0) + (l) + (2) 



MI2L1 



(93) 



(94) 



where the weights are given by (p 



3. Configurations {Si,S2,M) and (S'a , Si , M) 

The tree traps are the main trap M, the trap defined as before as the biggest trap before M, and the second 
biggest trap before M, that we call 5*2 if its position is between 5*1 and M, and that we call S2 if its position is 
between and 6*1. 

For the configurations (5*1, S2, M), the joint distribution of the rescaled positions and trapping times is given by 



Ts,^S2,m{XuX2,X;ts„ts,) = e{Ts, < rg, < t)e{0 < X^ < X2 < X) 

-xi 



TSi \TSiJ TS2 \TS2J 

and the corresponding diffusion front reads 

-PsistM^^^^(") =Pl(*;'rSi)^m,nSi + P2 (i; , Ts JJ™,„s^ + P3(i; Tfii , rsj(5™,„j,, (96) 

where the weights are given by (^7|). 

For the configurations (6*2, Si, M) the joint distribution of the rescaled positions and trapping times read 

Ts'^,Sum{X2,Xi,X;ts„ts„) = 0(ts, < ts, < i)0(O < X2 < X^ < X) 



fi f t y -X 



TSi \TSi / TS2 \TS2 

and the corresponding diffusion front reads 



(97) 



where the weights are given by (p 



4-. Configurations (Si , AI, Li 



The th ree traps are the trap 5*1 introduced in (QIC), the main trap M and the next renormahzed trap Li introduced 
in (IIIB). The joint distribution of the rescaled positions and trapping times is given by 



TsuM,L^{Xs,X, Xl;ts, Tm) = e{TM >t>Ts> l)e{XL >X>Xs>0) 



1^ f t 



The corresponding diffusion front reads 

where the weights are given by (^7|). 

We now use the statistical properties of these three-traps configurations to compute observables at order /i^ . 



(99) 



(100) 
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B. Diffusion front at order 

1. Contributions at order jj.^ of the two-traps configurations 

We have already studied the contributions of two traps configurations when studying the order /i. The contribution 
of order /x^ of the configurations MLi (|55|) reads 



Jo ^ 

Similarly, the contribution of order of the configurations SiM ( ^9|) reads 



V 



(101) 



(102) 



2. Contributions at order fi^ of the three-traps configurations 

The specific contribution at order of the three-traps configurations of type ML1L2 can be obtained by subtracting 
from (|9^) the two-traps configurations MLi (Esh 



j(0) + (l) + (2). 



,(0) + (l) + (2) 



The average over the samples with the measure (91) yields the correction of the scaling function ( p36| 

5mlil2(^) ^ y y '^'^^ J '^'^'^ J J Tm,Lx,L2{X,Xmli,Xl2;tm,tl^) 
P3{t; TM, TL,) id{Y - X2) - S{Y - X,)) 

Jo V Jo w 



2 



(103) 

(104) 
(105) 
(106) 



Similarly, the specific contribution at order fi^ of the three-traps configurations of type MI2L1 can be obtained by 
subtracting from (^) the two-traps configurations MLi ( ^ ) and this yields after averaging over the samples with 
the measure (p3|) 



P, 



(2) 



,(0) + (l) + (2) 



MI2L 



The correction to the scaling function ( pq ) thus reads 



5m/2Li(-^) 















2 


/o 1' Jl 



+00 



(107) 



(108) 



For the three-traps configurations of type S1S2M , we have to subtract from ( p6| ) the two-traps configurations SiM 
(^l|) and to average over the samples with the measure (95) 



The correction to the scaling function ( p6| ) thus reads 



3(0) + (l) + (2), 



j(0) + (l) + (2). 



y2- 



/ -/ ^P2(x',U;)-fO(M^) 



(109) 



(110) 



For the three-traps configurations of type S^SxM^ we have to subtract from (p8| ) the two-traps configurations S\M 
( pT| ) and to average over the samples with the measure (97) 



Pl{t;TS2)Sm,,l.s' + iP2it;TS2,TS,) -p2(i;0,TsJ)5^,„s^ {p3it;TS2,TsJ - P3 (* ; 0, Ts J )5„i^„j 



(111) 

(112) 
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The correction to the scaling function (Bq) reads 



' — In we 

V 



Y2 r+oc , p+oc J 

e-^(y-— )mM -/ —e{v<w)[p2{v,w)]+Oif,'){U3) 
2 Ji V Ji w 



For the three-traps configurations of type SiMLi, we have to subtract from ( |100D the one trap configuration (|3C|), 
and the corrections due to the two-traps configurations MLi (52)and SiM (|5^), and to average over the samples with 
the measure ( p9|) 



H - P^^^^r^ in) - in) - 4:^'^^ in) - P^^j^!'' in) 

with the scaling function 



[pi{t; , tm) - Psit; 0, rM)](^^ 



The sum of all contributions of order /i^ finally reads 







dw r 








2 


ll w Jo 



-P2{v,w) +0{n^) 



(114) 



(115) 



gllUy) ^ iY) + gZiL. (Y) + gf^s^MiY) + Ss^^iY) + gf^niY) + 9Z{iL2iY) + 9^l. iY) (US) 



dw 



P2iv,w) 



+^i^e'^{Y^ - 3r + 1) [ — liive-^ - f — lnt>(l - e"") 
Ui ^ Jo ^ 

The double integral may be computed as follows 



dw 



'-+°° dv /■+°° 
and we obtain the final result 



P2{v,w) 

W ./q w 



dw 



dz- 



1 



z- 1 



dz 





(Inz) 



TT 



(117) 



(118) 



2 12 / I 2 12 



y2 f ll _ ^ 
2 12 



which coincides with the expansion (nil) of the exact diffusion front described in Appendix 



(119) 



C. Thermal width at order ^ 

1. Contributions at order /i^ of the two-traps configurations 

We have already studied the contributions of two traps configurations when studying the order ^. The contribution 
of order /i^ of the configurations MLi reads (p^) 



Jo ^ 

whereas the contribution of the configurations SiM reads ( |65|) 

[A(M)]i''M = -6M' / -ln«e-"(l-e-'') 



(120) 



(121) 
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2. Contributions at order of the three-traps configurations 

For a given configuration of three traps situated at (ni, rj2, ^3) with occupation probabihties (pi,P2,P3), the thermal 
width reads 

< An^(i) _PlP2(f^2 - nif + pipsin^ - nif + P2Pz{n3, - 712)^ (122) 

Following the procedure described above for the diffusion front, we obtain the specific contributions at order /i^ of 
the various configurations as follows. 

The configurations of type ML1L2 with the measure ( pl| ) give 

[ A(A^)]i^A,^, ^ [A(A^)]iyz!'/f! - (123) 

= P2[t\TM ,TL^)pz{t;TM ,tli)[Xl2 -Xl^Y + Pi{t\TM)pz[t\TM ,tl^)[[Xl^ - XmY - (Xli - Xm)^] (124) 

= ^2 ^ —[2p2{v,w)p3{v,w)+4pi{v)p3{v,w)]+0{fi^) (125) 



The configurations of type M/2L1 with the measure ( p^ ) give 

[A (M)]i^U: ^ [M^^)]^n.L^. " [Mp)]i]^L!'^ (126) 

= P2{t;TM,Ti^)p3{t;TM,Ti2)[XLi - Xi^Y +pi{t\TM)P2{t\TM,Ti^)[{Xi^ - XmY - {Xl^ ~ XmY] (127) 

= /x^/ — / —[2p2{v,w)p3{v,w) ~ Api{v)p2{v,w)] (128) 
Jo V Ji w 

The configurations of type S1S2M with the measure (^) give 

[ A(M)]i%,M ^ - [a(m)]1,?;/'^ (129) 

= P2{t; Ts,,TsMt; ts„ts,)[Xm - Xg^]^ +pi{t; rsMt; ts,,ts,)[{Xs, ~ Xs.Y - {Xm - Xs,Y\ (130) 

= A* / — / — [2p2(w,ti')P3(w, w) - 4pi(w)p2(w,ti')] (131) 



w 

The configurations of type S'^S\M with the measure ( |97| ) give 

l(2) - rA/„M(l) + (2) rA/„M(l) + (2) 



[A(m)]^Wm ^ [A(m)]^7;m^ - [A(a^)]^Vm (132) 
= Pi (i; )vi {t; TS2 , ) {Xs, - Xgry + pi {t; , )P3 {t; ts^ , ts^ ) {Xm ~ ^s; )^ 



+[P2{t; , )p3(tTs2 ,TsJ - P2it;0,Ts^)p3{t;0, Tg^ )] (^Af - (133) 
= ^2/ — / — [2pi{v)p2{v,w) +6pi{v)p3{v,w) + 2p2{v,w)p3{v,w) ~ 2pi{w){l - pi{w))] (134) 



1 ^ ^1 



The configurations of type SiMLi with the measure ( p9| ) give 

[ A(A^)]i%/L, ^ [A(M)]i^?;A',^ ~ [A(A^)]i\^;/^^ - [Mp)]^^^^ !'^ (135) 
= PI {t; Ts, )P3 {t; Ts, , tmWXl, - Xs,y - [Xm - J^] 



+ [P2it;Ts^,TM)p3{t;Ts^,TM) ~ P2{t;0,TM)p3{t;0,TM)]iXL^ - XmY (136) 

= — / — [2p2{v,w)p3{v,w)+Api{v)p3{v,w)~2pi{w)[l-pi[w))\ (137) 



Finally, the sum of all contributions at order reads 

[A(M)]il, ^ [A(M)]g5.M + [A(A^)]i!^M + [A(m)]L%,m 

+ [A(M)]i''ML, + [A(m)]£'z.l, + [A(A^)]i^'i,i, + [A(A.)]i^)^^ 



Jo 

p-\-oo p-\-oc 

— —[^p2{v,w)p3{v,w)-Qp^{v)p2{v,w)+2pl{w)p3{v,w)]+0{^J!^) (139) 
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The double integral may be computed as in (118) and yields 



— / — [4:p2{v,w)p3{v,w) ~ 6pi{v)p2{v,w) + 2pi{w)p3{v,w)] ^ 41n2 (140) 

V ./„ w 6 



and thus the final result 



coincides with the expansion of the exact result (ha 



[A(m)]JoL = M'[21n2(ln2 + 2^e) - - - 41n2] (141) 



D. Localization parameters at order fi^ 



1. Contributions at order of the two-traps configurations 



We have already studied the contributions of two traps configurations when studying the order /i. The contribution 
of order /i^ of the configurations MLi (^) reads 

[Yu]%, =f^' f'-lnv [e-'^^ + (1 - l] (142) 

Jo ^ 

whereas the contribution of order fj,^ of the configurations SiM ( ]69| ) reads 

[i^fcls.M = / - invie-"^ + (1 - e-n" ~ 1] (143) 

2. Contributions at order jj.^ of the three-traps configurations 



For a given configuration of three traps with occupation probabilities (|87|) , the localization parameters read in terms 



of the variables v = — and w = — 



Yk=p'^iv)-hP2{v,w)+p^{v,w) (144) 



Following the procedure described above for the diffusion front, we obtain the specific contributions at order of 
the various configurations as follows. 

The configurations of type ML1L2 with the measure (^ij) give 



9 dv dw 



[pI{v, w) + pI{v, w)-{1- p^{v)f] + 0{i?) (145) 



The configurations of type MLI2L1 with the measure (93) give 



in&L. - - in]^^^r^^' (146) 

= m'/-/ —[p'2{v,w)+pl{v,w)-{l-p^iv))'']+0{f,') (147) 
Jo V Jl w 

The configurations of type S1S2M with the measure (p)5|) give 



rv;i(2) _ rv,i(")+(i)+(2) „ rv,i(o)+(i)+(2) 

= m'/ -/ — [p2'(«,^i')+P3(^,^)-(l-PiM)1+0(/i^) (148) 
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The configurations of type S^^SiM with the measure (^) give 



W _ \V l(0) + (l) + (2) _ \Y l(0) + (l) + (2) 



The configurations of type SiMLi with the measure (|9T|) give 



= ^^' - —[pl{v,w)+pl{v,w) + l-{l-p,{v)f-PiH~{l~Pi{w))''] + 0{fi') (151) 

Ji y Jo w 

The sum of all contributions of order fi^ thus reads 



= m'/ -/ — b^(z;,u;)+P2^Ki;)+2p^-KH + l-P^(^')-2(l-pi(«))'=-(l-pi(«;))1 (152) 



For the special case fc = 2, we find 



4 In 2 

6 



in agreement with the expansion of the exact result (p^). For the special case fc = 4, we find 



(153) 



[yitli = ( - Y + 2(21 ln2 + In3(ln3 - 12)) ) (154) 

E. Correlation function at order fi^ 

1. Contributions at order /i^ of the two-traps configurations 

We have already studied the contributions of two traps configurations when studying the order /i. The contribution 
of order of the configurations MLi (fz^ ) reads 

Pm(A)]^A, - /' - ln«e-'^(l - e-") (155) 

Jo 

whereas the contribution of order of the configurations SiM (|8l|) reads 



[C/.(A)]g^M - - VAe"^ / ^ln7;e-''(l-e-) (156) 



2. Contributions at order fi^ of the three-traps configurations 

For a given configuration of three traps situated at (ni,n2,n^) with occupation probabilities {pi,p2,P3), the two- 
particles correlation function reads 

C{l,t) = (pI +pI +pl)5ifi + 2pip25i 

,n2 — "1 

+ 2P2P3S1 (157) 

Since the weight of the delta pic is given by the localization parameter Y2 that we have already considered above, 
we will consider in the following only the scaling function C^(A) (^8|). Following the procedure described above for 
the diffusion front, we will obtain the specific contributions at order fi^. 
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The configurations of type ML1L2 give the specific contribution at order /x^ 

= e-^M^ f f dw p^^^^^ w)p3{v, w) - 2pi{v)p3{v, w) + X2pi{v)p3{v, w)] + 0{iX^) (158) 
Jo V Jq w 

The configurations of type MI2L1 give the specific contribution at order /x^ 

-A 



e 



— —[2pi{v)p2{v,w)+2p2{v,w)p3{v,w)-2pi{v)p2{v,w)\\+0{^i^) (159) 

Jo V Ji w 



The configurations of type S1S2M give the specific contribution at order /x^ 



= e- V 



/ — / — [2pi{v)p2{v,w) +2p2{v,w)p3{v,w) -2pi{v)p2{v,w)X]+0{fj,^) (160) 
Jl V J^ w 



The configurations of type 6*2 5*1 M give the specific contribution at order fi"^ 

'• + OC J,, /•+00 



= e-^/i^ / — / — [2pi{w)p2{w, v) + 2p2{w, v)p3{v, w) - 2pi{v){l - pi{v)) + 2pi{w)p3{v, w)X] + 0((/tl%)L) 

Jl V J„ w 



2 



The configurations of type SiMLi give the specific contribution at order /u' 

= H^e-^ — — [2p2(w,t^)P3(t^,w)-2pi(w)p3(w',v)-2pi(w)(l-pi(t;)) + 2pi(w)p3(t;,u;)A]+0(A/(^)62) 

Jo Jl ^ 

The sum of all contributions of order fi^ reads 

M^e-^ y-ln2(4 + ln2 + 7B) + A(-y +In2(ln2 + 7E)) (163) 



V. HIERARCHICAL STRUCTURE OF THE IMPORTANT TRAPS 

It is now clear that the procedure that we have described up to order /i^ can be generaUzed at an arbitrary order n 
: all observables at order /x" can be obtained by considering a dispersion of the thermal packet over at most (1 + n) 
traps, that have to be chosen among a certain number f2„ of possible configurations of the traps. Our aim in this 
Section is not to pursue any further explicit computations but to get some insight into the set of important traps that 
play a role at a given order n. 

A. Set of the important traps at order n 

At order n, the important traps are : 

• the main trap M 

• the following n large renormalized traps Z/i,...L„ 
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• the n biggest traps Si...Sn among the small traps before M 

• the {n — 1) biggest traps I2'^\..I^^ among the small traps in the interval between M and Li 

f21 (2) 

• the (n — 2) biggest traps /g ...In among the small traps in the interval between Li and L2 

• ... 

• the biggest trap /i" among the small traps in the interval between L„_2 and £„_i. 

The index at the bottom represents the order of occupation in fi as in the Figure |^. The total number of traps is 
thus 

T^ = l + n + ±^ = l + ^^^ (164) 

1=1 

which generalizes Ti = 3 (M,Si,Li) and T2 — 6 {AI, 81,82 identified with 82,Li,L2)- 

B. Set of the important configurations at order n 

With these r„ traps, we have now to construct the possible ri„ configurations of (1 + n) traps, that are ordered by 
in positions, and that contribute up to order /u". We have 

n 

r2„ = rj„-i + cj„ = ^tjj (165) 

i=0 

where w„ represents the number of configurations that begin to contribute at order n. 
We may now decompose 



= a(^")+a(^"-) + ...a^)+af (166) 
lat coi 

configuration M LiL2...Ln, whereas for j = 



where a^^^^ is the number of configurations that contain Lj as the rightmost trap. For j = n, there is only ai^"'' = 1 



= n\ (167) 

since we have to order in space the n traps 81... 8n before M. More generally, at order j, to construct the configurations 
of (n + 1) traps containing MLi...Lj, which represent (j + 1) fixed traps, we have to choose (n — j) traps among the 
(j + 1) available intervals and to count the possible positional orders in each interval 

+00 -\-oo 

«n^ = E- E S{Y,P^^n~J)pl\...p,+l\ (168) 

pi=0 Pj + i=0 i = l 



The final result is thus that the number of new configurations that appear at order n reads 

(169) 



^" = E 

J=0 



+ OC +00 j+1 

pi=0 pj + i=0 j=l 



which generalize what we have found before for the lowest orders ujq — 1 (M), uji — 2 {81M and MLi) and uj2 — 5 
{8182M, 8'28iM, MIiLi, ML1L2 and 81ML1). 
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VI. QUANTITATIVE MAPPING BETWEEN THE BIASED SINAI DIFFUSION AND THE 

DIRECTED TRAP MODEL 

A. Renormalized landscape for the biased Brownian motion 

The real space renormalization group (RSRG) method can also be applied to the biased Brownian lands cap e | 
The distributions of the barriers F = F + ^ in the renormalized landscape at scale F are given by |l^, 



1 

_ 26 
1 ~ 



PfiO =^—^e-^TZ7-^^2Se-^'' (171) 



where the parameter 25 reads in terms of the notations (|^, ||) 

2^-^ = ^ (172) 

1 (7 

As a consequence, the distribution Py {C) of barriers against the bias can be considered as infinitely large only in 
the limit of vanishing bias J — > 0. It is only in this limit that all particles of the same thermal packet remain in the 
same renormalized valley asymptotically. 



B. Trapping-time of a renormalized valley of barrier F 

Let us now recall a standard result for one-dimensional Fokker-Planck equation ||2l[| : for a particle diffusing in a 
potential U{x) on a interval [a, 6] with reflecting condition at a, the exit time defined as the first-passage time 0{x) 
at the point b for a particle starting at x € (a, b) at time t = Q can be studied for an arbitrary potential U {x) : the 
moments 



are given by the recurrence 



%{x) =< [e{x)T > (173) 



9,,[x) = fi f dye^'^^y^ r dze-^^^'-^nen-iiz)] (174) 



with the initial condition 9n=o{x) = 1. In particular, the first moment reads 

ei{x)^(3 f\ye'^^^y'> r dze~^^^^'-^ (175) 

For the biased Brownian landscape, the exit time over a barrier F when starting at the bottom of a renormalized 
valley that we choose as the origin can be obtained (see Figure ||) by choosing a at the height F on the renormalized 
descending bond on the left and 6 at a potential F after the top of the barrier F. It seems that usually [|l] one chooses 
b exactly at the top of the barrier to derive the Arrhenius factor, but we think that to obtain the correct prefactor, 
one has to choose b on the descending potential after the top to be sure that the particle will not return in the trap 
where it started. Indeed, when the particle sits just on the top, there is a finite probability to return to its starting 
trap, which is for instance a probability 1/2 for a potential that is symmetric around its top. So for a given realization 
y of a renormalized valley, the first moment of the escape time reads 



b 



9,{V} = p / dyeP^'^y^ 



i-y 



(176) 

where the biased Brownian potential V{x) = —Fqx + U{x) satisfies the constraints of a renormalized valley at scale 
F (see Figure |^) : it starts at 1^(0) = e, it then evolves on each side x > and a; < in the presence of absorbing 
boundaries at and F, and is conditioned to finish aX V = T and not a,t V = 0. On the negative side, {x = a) is the 
random position where the potential first hit F. On the positive side, after the random position Ir where the potential 
first hit F, the potential again evolve in the presence of absorbing boundaries at and F, and is conditioned to finish 
a.t V = 0, at some random position called b, and not at F = F. 
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FIG. 2: Computation of the escape time fro m a renormalized valley of barrier F : we consider the first-passage-time at b for 
a particle starting at 0. The double integral (L75) is dominated by the Arrhenius factor e^^ , and the prefactor is the product 
of two partition function : Zv represents the partition function of the bottom of the valley and Zb represents the partition 
function of the inverse potential ( — V^) near the top of the barrier T. 



^P^'' where 



As usual for the problem of escape over a large barrier, the double integral ( |176| ) is dominated by the saddle where 
V{y) is maximal and where V{z) is minimal. For a renormalized valley (see Figure ||), these regions are y 
V(y) ^ r and z where V{z) ^ 0. This saddle-point analysis yields the following leading behavior 



0^{V} 

The prefactor is simply given by the product 

to{V)^PZbZv 

where Zy is the partition function of the infinitely deep renormalized valley 



Zv = lim 

r^oo 




where the random potentials 



V+ix) 



-Fox + Ui{x) 
Fax + U2ix) 



(177) 



(178) 



(179) 



(180) 
(181) 



are defined in terms of two independent Brownian trajectories Ui{x) and U2{x) (||) starting at V+{0) = e = V-{0). 
The potentials V±{x) evolves in the presence of absorbing boundaries at and F, and are conditioned to finish at 



V = T and not at = 0. ^p^^ and V^' are the random times where V± respectively first hit V = T. 

Similarly, the factor Zg is the partition function of an independent infinitely deep renormalized valley, which 
represents what happens in the vicinity of the top of the barrier F when considered with the change V —V to 
transform it in a valley (see Figure]^). For the biased Brownian landscape considered here, by symmetry, Zb is simply 
an independent realization of the variable Zy- 



,(2) 



The same saddle point analysis may be applied to higher moments given by the recurrence (174) to obtain 



(182) 



So for a given renormalized valley of barrier F, the escape time t is distributed exponentially as in the trap model 

1 



/ei{y}(0 



e,{v} 



(183) 
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where the trapping-time 6'i{V^} (177) depends mostly on the barrier F via the usual Arrheniu s fac tor e^^ , but also on 
the details of the structure of the valley near the bottom and near the top via the prefactor (178). 



C. Distribution of the trapping-time of renormalized valleys 



We are now interested into the distribution of the trapping-time 6*1 {V^} over the ensemble of reno rmalized valleys 
existing in the renormalized landsca pe a t scale T. The distribution of the barriers is given by (171). So we have to 
study the statistics of the prefactor (178). 

It is more convenient to work with dimensionless quantities by rewriting the partition functions as Zy = zi/{(t0^), 
Zb = Z2l{p0^) so that 



ro(^) 



(184) 



where z\ and Z2 are independent random variables whose probability distribution P{z) is characterized in Appendix 
P by its Laplace transform 



/.-l-oo 

/ dze-'''-P{z) 



1 



r(i + M)/p(2^/i) 



(185) 



where is the Bessel function of index /i. In the renormalized landscape at scale F, the probability distribution of 
the trapping time r of the renormalized valleys thus reads using (171) 



+ 00 /' + 00 



+ 00 



dziP(zi) / dz2P{z2) S 

"'0 "'0 

dz^z'^Pizi) / dz2iz2rPiz2) 



p/3(r+?) 



-Z1Z2 



(186) 
(187) 



So we have to compute the non- integer moment of order /i of the variable z. Using the integral representation valid 
for < ^ < 1 



+ 00 



r(i - /i) Jo 

we obtain the moment from the Laplace transform (|185|) 



dss-i-^ (1 



+ CX) 



dzz^P(z) = 



+00 



r(i - lA Jo 

r(i-M) 



1 - 



1 



dw 



w 



T{1 + fi) I^{2^s) 
1 



F2(l -|-^)4A'w/2(w) 



Using the wronskian property of Bessel functions, and their series expansion at small argument, we finally get 



+00 



dzz^'P{z) 



1 



s lim 
F(l - fl) a^O 



2^ K^{a) 



1 



r(i + M) 



(188) 

(189) 
(190) 

(191) 



The final result is thus that the distribution of trapping time r of the renormalized valleys existing at scale F reads 



Prir) =^-' ' 



/3r 



1 



T\a'^P^Tj F2(l-h/Li) 



(192) 



D. Precise choice of the renormalization scale F as a function of time 



We have seen in the trap model that the distribution of renormalized traps at t reads (^ 



qt{T) = eit < r)^ - 

T \ T 



(193) 
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To make it exactly coincide with the distribution (192) of the biased Sinai model, we have to choose the renormalized 
scale r of the landscape to be the following function of time 



r(t) = Tin ta^P^ (r2(l + ^ 



(194) 



The RSRG method gives that the distribution of the length 1+ of the descending renormalized bonds is simply 
exponential in the limit F — )■ oo 



(195) 



where the mean length reads 



25r 



2Sr 



so that it reads as a function of time (194) 



(196) 



(197) 



L ^' J 

This length scale b{t) exactly corresponds to the ratio ^ '''°''-^(^^) of the constants appearing in the exact diffusion 
front of the two models (A26). 



E. Usual RSRG in the limit -> 

It has been shown in ||^ that the "effective dynamics" , where at time t, the particle is typically at time t around 
the minimum of the renormalized valley containing the initial condition, is sufficient in the double limit t ^ 00 S ^ 
with 



7 = Sr{t) (198) 



fixed and X = fixed. 



The limit 7 — > corresponds to the symmetric Sinai diffusion, whereas in the limit 7 — > 00, the model becomes 
directed at large scale and the diffusion front converges towards 12] 



P{x,t\Q,Q) ~ e{x)—e-Wl (199) 
7^00 bit) 



where h{t) represents the mean length of renormalized descending bonds ( |197| ) . This limit actually corresponds to the 
limit /i ^ of the exact Levy front jlj, |l^, |l^ as described in the Appendix (^). 



Spreading of the thermal packet over many renormalized valleys 



The renormalized valleys of the Sinai model with bias are the analog of the traps in the directed model . For /i — > 0, 
the bottom of the renormalized valley containing the origin described above (199) is the analog of the main trap M 
described in section ||. 



At first order in /i, as in Section III A, there are two effects : 



• Next renormalized valley Li 
The main renormalized valley M at scale V{t) has a tra ppin g time tm which is distributed as in the directed 



trap model (|40|), since we have defined the relation T{t) (194) by identifying the trapping time distributions of 

the two models. So there is a small probability (1 — e ) that the particle has already escaped from this main 
renormalized M at time to jump into the next renormalized valley Li. 
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FIG. 3: Hierarchical structure of the important valleys for a particle starting at the origin. The barriers against the bias that 
are emphasized by the straight lines correspond to the depths of the trap model represented on Figure ^. The bottom M of 
the renormalized valley that contains the origin at scale V is occupied with a weight of order 0{fi^). The bottom Li of the 
next renormalized valley and the bottom Si of the biggest sub- valley before M are occupied with weights of order 0{fi). The 
next-nearest renormalized valley L2, the biggest sub- valley I2 between M and Li, and the second biggest sub- valley S2 before 
M are occupied with weights of order 0{p?). 



Moreover, the RSRG approach |^ yields that the joint distribution of the trapping time tm and of the positions 
xm and xl of the bottoms of the main renormalized valley and of the the next renormalized valley Li reads 

Vm,lAx,xl;tm) =0{t<TM)OiO<x<XL)^Tj7-e (200) 

= eit<TM)e{0<x<XL)^T^e-WJ (201) 



which is the analog of (|39|). The only change is in the prefactor in front of i'' in the scale b(t) (197). 

• Last decimated renormalized valley Si 

The last decimated barrier against the bias inside the main renormalized valley between the origin and the 

bottom, defines the last decimated sub-valley Si : it has a trapping time ts^ < t which is not zero and thus 

_ t 

there is a small probability e that the particle is still trapped in the sub-valley Si at time t. 

Moreover, the RSRG approach yields that the joint distribution of the trapping time ts and of the positions 
xs and X of the bottoms of the last decimated valley 5*1 and of the the main renormalized valley M reads 

'Ds,.m{xs,x;ts) ^0{t> Ts)e{0 <xs <x)^—^—e-^ (202) 

TS b[t)b[Ts) 



which is the analog of (|45|). The only change is again in the prefactor in the scale b{t) ( |197| ). 
It is clear that this analysis may be generalized to further orders in /i. 



G. Conclusion : Equivalence of the two large-scale renormalized descriptions 

The statistical properties of the spreading of the thermal packet over many renormalized valleys and sub-valleys 
inside the main one are thus exactly the same as in the directed trap model discussed in details in the first Sections. 
In particular, the localization parameters Yk of the trap model represent coarse-grained localization parameters for 
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the biased Sinai diffusion : "at the same position" in the trap model means "at a finite distance around the bottom 
of the same renormahzed valley" for the biased Sinai diffusion. As a consequence, for all rescaled quantities 
the results are exactly the same up to the global prefactor in the scale b{t) (197) : this was already known for the 
averaged diffusion front (see Appendix^), but this also holds for the thermal width (^7|) , for all other rescaled thermal 
cumulants (24), and for the long-range part of the two-point correlation function (^9[). 

The new property of the Sinai model is thus the internal structure of a renormahzed valley, that induces a dispersion 
over finite distances of the particles that are in the same renormahzed valley. We now study the statistical properties 
of the biased Brownian valleys. 



VII. INTERNAL STRUCTURE OF THE TRAPS IN THE BIASED SINAI DIFFUSION 
A. Probability distribution inside a renormalized valley 

The probability distribution of particles inside the same renormalized valley can be obtained by generalizing the 
approach of the Sinai symmetric case Q : for each realization of a renormalized valley, it is given by the Boltzmann 
distribution on this valley. So asymptotically as i ^ oo, the probability distribution of the distance y to the bottom 
of the valley, averaged over the environment reads 



Pviy>o) 



lim 



-0V+iy) 



{V+},{V-} 



Py{y < 0) = lim 

r^oo 



^-f3V^{\y\) \ 



(203) 



where the random potentials V± satisfy the same conditions as in ( |l8l|) . 

The computation of the functional (202) is given in Appendix M. It yields the non-intuitive result that the 
probability distribution Pviy) is actually symmetric in y — > —y. The restoration of this symmetry comes from the 
conditioning of the biased random walk to reach T on each side. Its Laplace transform r eads ( pl^ ) 



Pv{p) = 



due-PyPv{y) 



1 



ds- 



,0 r^i + ti)J, Us)i,is)J, 

where the only factor containing the Laplace parameter is the index 

2T2p 



dzzl^{z) 



(204) 



(205) 



(206) 



cr cr/i 

For fi > 0, the series expansion in the Laplace parameter p is thus regular leading to 

1 2T-^?7 
Pvip^O) ^D{^,) + ... 

2 fia 

All moments are thus finite, contrary to the symmetric case /z = |5| , where the behavior as (1 — y/pc+ ..) corresponds 
to the algebraic decay as Xjy^l"^ . So in the biased case, the distribution inside a renormalized valley is very narrow, 
contrary to the symmetric case. 

B. Localization parameters inside a renormalized valley 



For fc particles that are in the same renormalized valley, the localization parameters may be computed as an average 
of the k-th power of the local Boltzmann weight over the infinitely deep biased Brownian valleys ( |181[ ) . Generalizing 
the approach of ||^ to the biased case, we have 



£=± 



-t-oo 



y— , 



r' + OO . . / p-\-00 \ 



(207) 



(208) 
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Using the results ( |B18 ) of the Appendix, we finally get 



.(§) 



2(1-1 



IM 



(209) 



C. Correlation function inside a renormalized valley 



The correlation function of two particles at Boltzmann equilibrium in an infinitely deep biased Brownian valley 
reads 



Cvalleyil > 0) = 2^ / dy 



-pvAy)-0v,(y+i) 



^-0V+(v)-l3V_(l-y) \ 



-2 dy( 



where the average < .. > is over the realizations (V+,y_) satisfying ( |18lD 
Using the explicit results of Appendix O, we finally get 



Cvalleyiv) — p2Q 



ds (f) 



2fi 
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ds 



\ 2m 



T^l + ^,)Jo s IM 



dzziM[KM-^^iM 



(210) 



(211) 



VIII. DISCUSSION OF THE UNIVERSALITY 



We now briefly discuss the question of the universality. The RSRG method that describes the large scale structure 
at scale F of the random potential is valid for all discrete models wi th r andom forces The parameter 25 that 
describes the distribution of barriers against the drift at large scale (171) may be expressed for a discrete random 
force model as the non-zero solution of the equation p| 



-2.5/ ^ 1 



(212) 



which is known to determine the anomalous diffusion exponent /i — 26T [ |15| , |lq | . So for a given value of the parameters 
(26, cr), the renormalized landscape at scale T is universal. 



However, it is clear from the analysis of the escape time of a renormalized valley (177) that the prefactor ( 178 ) in 
front of the Arrhenius factor e^^ is not universal : the partition functions Zy and Zb depend on the details over 
finite scales of the potential near a bottom of a renormalized valley and near a top of a barrier. 

So for a potential that belongs to the universality class (25, cr), but that is not a biased Brownian at small scales, 
the distribution of the trapping times in the renormalized landscape at scale F reads 



Prir) 



^ (z(;)(z^) 



(213) 



so that the quantitative mapping onto the trap model (193) is realized for the choice of the RG scale F as a function 
of t according to 



T{t) = Tin 



which corresponds to the length scale 



b[t) 



(214) 



(215) 
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This shows that the factor //^ is universal and comes from the mean length of descending bonds in the renormalized 
landscape at large scale, whereas the factor r^(l of the biased Brownian motion is not universal and comes from 



the probability distribution of the partition function of a biased Brownian valley (191). However it is expected to 
be valid for discrete models in the limit where the lattice constant is very small as compared to the thermal length 
It = T'^/d- For the localization parameters and the correlation function of two particles inside the same renormalized 
valley, the discussion of the universality is the same as in the symmetric case [g . 

IX. CONCLUSIONS AND PERSPECTIVES 

To study the anomalous diffusion phase a; ~ i'^ of the directed trap model and of the Sinai diffusion with bias, we 
have extended the usual RSRG method that assumes a full localization in a single valley to allow for the spreading of 
the thermal packet over many renormalized valleys. We have shown how all observables can be computed via a series 
expansion in /i : at any given order /x", it is sufficient to consider the spreading over at most (1 + n) traps. We have 
given explicit rules for the statistical properties of these traps. We have shown the exactness of these expansions in 
/i by comparing up to order n = 2 with the already known exact results, such as the diffusion front [l^ , the thermal 
width [|l9j and the localization parameter Y2 j2^ . Our construction moreover gives a clear physical picture of the 
localizations properties in the anomalous diffusion phase, and explains the typical shape of the diffusion front in a 
given sample obtained by numerical simulation (Figure 4 of pO|). 

In a forthcoming paper ||2^ , we will adapt our method to study the localization properties and the aging behaviors 
in the symmetric (i.e. undirected) trap model which has attracted a lot of interest recently | p3| , p^ , p5[ |. 

For the field of biased diffusion in one-dimensional random potentials, it would be very interesting to study the 
influence of correlations on the localization properties studied here for the Brownian case. In particular, the case of 
algebraic correlations {U{x) — U{y)Y ^ \x — is known to give rise to a creep motion for < 7 < 1 [Q. For DNA 
sequences, it seems that the interesting cases are not only the Brownian case 7 = 1 |^ but also the values 7 > 1 [ p7[ . 
Another physically interesting case concerns the logarithmic correlations, which give rise to a freezing transition in 
the dynamics |Q as well as in the statics p9| . 

From the point of view of the RSRG method, since the usual RSRG is asymptotically exact for infinite-disorder 
fixed points the extension introduced here can be seen as a systematic expansion in the inverse disorder strength. 
It can therefore be used in the field of random quantum spin chains Jl3t to study the Griffith phases, and for the 
classical random field Ising chain in the presence of a small magnetic external field [Q . 

Finally, the expansion in the important traps for the dynamical models discussed in this paper has a static coun- 
terpart, with the following differences : in the static case, the expansion parameter is the temperature T, and the 
main trap M corresponds to the absolute minimum of the random potential. We have already shown in for the 
toy model consisting of a Brownian potential plus a quadratic potential, how the thermal cumulants at first order 
in T can be explained by studying the statistical properties of the configurations presenting two nearly degenerate 
minima. We will discuss in [ pi] in a more general context the structure of the low-temperature series expansions in 
some disordered systems. 

Acknowledgments 

I thank Pierre Le Doussal for helpful discussions at the early stage of this work and for his remarks on the manuscript. 

APPENDIX A: USEFUL PROPERTIES OF THE LEVY DIFFUSION FRONT FOR < a* < 1 

In this appendix, we recall some useful properties of the Levy distributions (jl^, |3^ and references therein). 

1. Definition and properties of one-sided Levy stable laws 

The rescaled sum 

1 " 

y = —Y.^r (Al) 
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of n identical independent positive random variables distributed with a law presenting the algebraic decay 

where < /i < 1, has for limit distribution as n — *■ cx3 the one-sided Levy law L^c(^L:A){y) defined by its Laplace 
transform 

r + <x> 

/ dye-^yi^Av) = e-^'" (A3) 

^0 



and where the constant c reads 



A) ^ (A4) 

' sm7r^r(l+/i) ^ ' 

In this paper, we will only use the following series representation ^ 

i^^Av) = --E (-4)' TTi^^--'^^ 
^yjr[\ J r(i + fc) 

which is convergent in the whole phase < < 1. 

We stress here that we have defined the constant c by the Laplace transform ([A^). Writing the inverse Laplace 
transform as a Fourier integral yields 



' 2nr ,1^ 27r 



so that the constant C appearing in the usual Fourier transform of Levy distributions 



, 27r 



reads in terms of the Laplace constant c 



C^ccos^ (A8) 



2. Levy diffusion front for the trap model 

For a given trap r, the distribution of the escape time t is exponential 



Mt) = -e-i (A9) 
r 



which yields after averaging over r 



TAt) ^ I dMr)Mt) = r ( -) . - ^^^^ (AlO) 



^0 



For a given sample (tq, ti, ...), the probability Pt{n) for the particle to be in the trap n at time t reads 

/+C>D 
Y[dUfrAU)d{to + tl... + t„^l <t<to + ti + ... + tn) (All) 

The average over the disorder 

Pt{n) = f f[dt~U{U)e{ti+t2 + ... + t„-i <t <ti+t2 + ... + tn) (A12) 



, +00 
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shows that the diffusion front at large time is dire ctly related to the properties of the sum of a large number n of 
independent variables distributed with the law ( AlO|) presenting an algeb raic decay ( [A^ ) : the rescaled variable 
y = ;ji77r is distributed with a one-sided stable Levy distribution L^^ct^apip.) {^)^ where the constant ctrapin) reads 
for the case (|A10| , |A^) 



TTfJ, 
sin TTfl 



The variable X = jj^ — y ''is thus distributed with the law 



1 



(A13) 



(A14) 



with the special value c = ctrapilJ^)- 



In particular, the series expansion (A5) gives the following series representation for the diffusion front 

/-^(^) = - L (-c^) YiiTW ^ (fc-i)!r(i-M 



Using the series expansion 



(A15) 



+ 00 



with the first coefficient 



r(i 



1 + 



771 — 1 



d2 



A. Ill 
2 12 



where je denotes the Euler's constant, we get the expansion in fi of the series (A15) for fixed c 

= ce-=^ [(1 - difi + d^fi^) + {difi - 3d2fi^)cX + d^ti^cX)^ + 0(m')] 
Expanding also in /x the value ( [aT^ ) 



(A16) 

(A17) 

(A18) 
(A19) 

(A20) 



ctrapif^) = 1 + -rM^ + O(m^) 
6 



we get that the diffusion front reads up to second order in fj. 



2 2 
g{X) = = [1 - diM + (rf2 + y )M']e"'' + MiM - (3^2 + ^)^i^]Xe~'' + d2/i'^'e-^ + 0(/i^) (A22) 



(A21) 



Using the numerical values (A19), we finally get the expression (|14|) of the text. 



3. Levy diffusion front for the biased Sinai model 



The exact form of the diffusion front was first determined in flStl for a corresponding discrete model. For the 
continuum model, the result has been proved in t he th eorem 1 of |T^, which states that, for < /i < 1, the rescaled 
variable X = has for probability distribution (A14), where the constant Cs{p) is given by a complicated implicit 
expression in . The value of this constant has been proven in p8[ to have the following simple expression 



TT/i 



2r2(/i) sin TT/i 



(A23) 



where we have used 
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In p4| , the same form was conjectured from the heuristic equivalence with the directed trap model via an indenti- 
fication of the parameters on some observable 

Cs^na^i^^) = — (^Y J,JT^. (A24) 

xi \ 2 / 1 ^(/i) sni nil 

where xi = 2T Vg a nd ri — Xi/{2Dq) in terms of the diffusion constant Dq in the pure case. This expression indeed 
coincides with ( |A23| ) for the units T = 1, a ~ ^ and Dq = \ used in |l8| . 
In the notations used in this article = 2^ this corresponds to 

CsinaiifJ-) = 7 o osNi/ "FT? ^ '■ (A25) 



To compare with the directed trap model, it is convenient to consider the ratio of the two constants (A13) 

Ctrap{^i) ~ (fT2/33)pr2(^) 



(A26) 



So beyond the natural dimensional factors, there is still a function r2(/i) between the two models, whose origin will 
be discussed in details in the text. 



APPENDIX B: STATISTICS OF THE INTERNAL STRUCTURE OF RENORMALIZED VALLEYS 

1. Distribution inside a renormalized valley 



To compute the functionals ( |203D , we generalize the approach developed in ||^ for the symmetric case /i = 0. We 
first exponentiate the denominator 



poo 

Pooiy>0) = dqR^{q)S+{y,q) 
JO 

Pooiy<0) = [ 

Jo 



(Bl) 



dqR+iq)S^{y,q) (B2) 

Jo 
where 

i?±(g)^(e-'/o^^--'''*'^')^^^ (B3) 
\ / {v±} 

These two functionals may be expressed as 

Rp{q) = AA± dl Imi (r - e, l\e) (B5) 

pT r + co 1 

Spiy, q) - M± du dl lun (r - e, Z|u)e-''"F± {u, y\e) (B6) 

in terms of the path-integrals 

fV(l)=u , 
JV(0)=UQ 

where the symbol 9[o.r]{V^(a;)} means that there are absorbing boundaries at F = and V = T. The expansion of 
the quadratic term of the measure yields 



j±j,](u,/|Mo) = e^^^'T^("-"°)F[o,r](u,Z|wo) (B8) 
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where 



JV(0)=uo 

(BIO) 

represents the analogous path-integral for the symmetric case. Its Laplace transform has been computed in Equation 
(B18) of g). We get 



(Bll) 



and the Laplace transforms with respect to y 



+00 



S^{p,q)^ / dye-PySpiy,q)^ 



in terms of the function 
2 



The normalizations A/± ^^re obtained with the conditions i?p {q 0) = 1 



A4 = 



(B12) 



(B13) 



(B14) 



We thus obtain that there is no dependence in the sign ± for the functionals Rp{q) and S^{p, q). As a consequence 
we get the non-intuitive result that the probability distribution Pv{y) is symmetric in y — > — y. The restoration of 
this symmetry comes from the conditioning of the random walk to reach F. We note that similarly, the distribution 
of the random times 1^ and are also the same, since we have, with the notations of j2j 



P]?(C = 0,s) _ Up{s) _ Vs + 5^ sinh 5T 



P± (( = 0, 0) (0) 5 sinh V^/ITP 

The Laplace transform of the distribution inside a valley thus reads 
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dzz- 



(B15) 

(B16) 
(B17) 

(B18) 



The final result is thus given by Equation (204) of the text 



2. Partition function of a renormalized valley 



We now consider the probability distribution of the partition function of a renormalized valley ( 17£ ) 

1 r+oo 



(B19) 



where the potentials satisfy the constraints (181) 
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Its Laplace transform can be directly expressed in terms of the functions (B3) which we have computed before 



dZvV{Zv)e-'i^^ = RUq)R^{q) 



r(i + /i) 



(B20) 



After the rescaling Zy = zi/{a0^), this corresponds to the result (185) given in the text. 
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